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1. Introduction and motivation 

Although several unquenched simulations of the QCD equation of state have appeared in re- 
cent years [p, lattice investigations of the thermodynamics of the Yang-Mills (YM) sector are still 
of interest from a fundamental point of view. In particular, the gluon sector is relevant for the limit 
in which the number of colors N tends to infinity, where analytical simplifications take place [0]. 
Furthermore, the large-Af limit is also a key-ingredient in the AdS/CFT conjecture ^ and in mod- 
els of QCD based on gauge-gravity dualities [Q]. These motivations led me to address a numerical 
lattice study of the equilibrium thermodynamics properties at finite temperature in SU(A') gauge 
theories [BJ|; similar works include refs. [Q]. 

In this contribution, we discuss some of the findings of ref. [Q], including, in particular, a com- 
parison with the improved holographic QCD (IHQCD) model recently proposed by Kiritsis and 
collaborators |^] (see also refs. ^ for related works). The IHQCD model is an AdS/QCD model 
based on a 5D Einstein-dilaton gravity theory with a particular Ansatz for the dilaton potential, 
which reproduces the leading terms of the SU(A') j8 -function at high energies, and linear confine- 
ment with a gapped, discrete glueball spectrum at low energies. The IHQCD model involves two 
free parameters, which in refs. [Q] were fitted to match the results of previous lattice calculations. 
In the context of gauge/gravity dualities, we also discuss the deficit of the entropy density s with 
respect to its Stefan-Boltzmann limit s , and a possible comparison with the AdS/CFT prediction 
for the large-Af limit of the ,jV = 4 supersymmetric YM theory []|] : 

- = | + ^C(3)(2A)- 3 / 2 + ... (i.i) 

So 4 jZ 

in a temperature regime where the Yang-Mills plasma, while still strongly interacting, approaches 
approximate scale invariance. 

Next, we investigate whether, in the temperature range studied, the trace anomaly of the deconfined 



SU(A') plasma A may receive contributions proportional to T 2 [10]: 

T* = T 2+B ' (L2) 

Finally, we present an extrapolation of our results for p / (N 2 T 4 ),A/ (N 2 T 4 ) , e / (A" 2 T 4 ) and s / (N 2 T 3 
to the N — > °° limit. 



2. Lattice setup 



We ran lattice simulations of SU(A') gauge theories with N = 3, 4, 5 6 and 8 colors on isotropic 
hypercubic lattices of sizes A^ 4 and A^ 3 x N t , with Af $ = 20 (16) for N = 3 (N > 3) and N, = 5. We 
used the standard, unimproved Wilson gauge action, and updated the configurations using heat-bath 
for SU(2) subgroups [11] and full-SU(A / ) overrelaxation AO). For SU(3), the physical scale was 



set using Sommer's parameter tq [13] while for SU(A / > 3) it was set by interpolating the string 
tension values taken from refs. [14], or (at the largest j3 -values only) using the method of ref. [|15|]. 
On the lattice, the trace anomaly A = e — 3p is proportional to the difference of the expectation 
values of the plaquette at T = and at finite T : 



dT T 4 



dloga 



((£/ n )o -(£/□), 



(2-1) 
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A / 7*, normalized 10 the SB limit of p / ~f 
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Figure 1: Comparison between the curves obtained using the improved holographic QCD model (yellow 
solid lines) and our lattice results for the trace anomaly, pressure, energy density and entropy density eval- 
uated in the various SU(N) gauge groups. All quantities are normalized to their Stefan-Boltzmann (SB) 
limits, except for A/T 4 , which is normalized to the SB limit of p/T 4 . 



The pressure was determined using the "integral method" [jig]: 



t-IL log s~ L log s 



a 4 N?N T Jfr 



(U n ) ), (2.2) 



according to the integration methods discussed in ref. [17]. Eq. (2.2) can be affected by finite- 
volume corrections [18] (especially for very high temperatures JT9| ] or small volumes [^]), but at 



the temperatures O.ST c < T < 3.4T C investigated here, these effects are screened [21] and negligible 
within the data precision jjZ"2]]. The energy (e) and entropy (s) densities were obtained as: e = 
A + 3p,s = (A + 4p)/T. 



3. Results 

Fig. [I] shows that the results for the rescaled trace, pressure, energy density and entropy density 
are very similar for all groups studied in this work 1 ; this makes it plausible that the QCD plasma 
could be described by models based on the large-Af limit. Fig. [j] also shows the comparison with 
the improved holographic QCD model discussed above: the solid lines denote the results obtained 
in ref. [0]; the agreement between our SU(A / ) results and the IHQCD model is very good. 

On the other hand, comparing the S\J(N) simulation results with predictions derived from the 
Jf = 4 SYM model, is less straightforward: the regime in which the QCD plasma is strongly 
coupled is far from conformality — see the left panel of fig. ^[ In our simulations the S\J(N) plasma 
approaches approximate scale invariance only at temperatures about 3T C ; at such temperatures, it 
is still far from the Stefan-Boltzmann limit (top right corner of the diagram), and still strongly 
interacting. Interestingly, in that regime the S\J(N) and the JV = 4 SYM entropy densities (both 
normalized to their values in the free limit) appear to be close to each other, if one uses the value 
of the renormalized 't Hooft coupling in the MS scheme as the A parameter of the supersymmetric 
model (right panel of fig. ||). Incidentally, we note that a comparison of JV = 4 SYM and full-QCD 



lattice results for the drag force also yields A ~ 5.5 [25]. Similar observations can be useful to pin 
down appropriate parameters for AdS/CFT models of the sQGP p6|]. 

Another issue we investigated is the possibility that the trace anomaly at temperatures close to 



T c may receive contributions proportional to T 2 of non-perturbative origin JIQ|] . The simulations 



This feature would be compatible with some quasiparticle model descriptions |E3h 
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Figure 2: Left panel: the equation of state, displayed as p(e), reveals the strong deviations of the SU(N) 
plasma from the dotted straight line corresponding to conformally invariant models, which is approached 
only at temperatures around 3T C ; the dashed line is the weak-coupling expansion for SU(3) taken from 
ref. [g4|. Right panel: the entropy density (normalized to its value in the free limit) as a function of the 
running 't Hooft coupling in the MS scheme; the dashed line is the corresponding prediction for strongly 
coupled JY = 4 SYM [g], identifying the A parameter with the renormalized YM coupling. 

seem to confirm that this may be a general feature of all the gauge groups studied in this work: the 
left-hand-side panel of fig. |3] shows that the results for A/T 4 may be compatible with the behavior 



described by eq. (|1.2|). However, we cannot rule out the possibility that the data may actually be 
described by a more complicated functional form (possibly involving logarithms of perturbative 
origin). Finally, the right-hand-side panel of fig. [3] shows an extrapolation of our results to the 



N — > °° limit, based on the parametrization for the trace anomaly given in eq. (CI) of ref. [28]. 

We plan to extend the present study of S\J(N) thermodynamics to other observables which 
could be compared with gauge/gravity predictions P9Q , and to renormalized Polyakov loops. On 
the other hand, to understand which non-perturbative features of the sQGP are directly related to 
the simplifications occurring in the large-N limit and which are not, it would also be interesting to 
address similar high-precision thermodynamics studies in models in lower dimensions or based on 



smaller gauge groups, for which powerful numerical algorithms are available [|30|]. 
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